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Abstract 

The averaged learning equation (ALEH) applicable to the prin- 
cipal component analyzer is studied from both quantum information 
geometry and dynamical system viewpoints. On the quantum infor- 
mation space (QIS), the space of regular density matrices endowed 
with the quantum SLD-Fisher metric, a gradient system is given as 
an extension of the ALEH; on the submanifold, consisting of the diag- 
onal matrices, of the QIS, the gradient flow coincides with the ALEH 
up to a local diffeomorphism. 

PACS: 02.40.Yy, 03.67-a, 02.40.Vh 

1 Introduction 

Quantum computing has been widely known to be a research area in rapid- 
rate progress. It has been recognized also as an interdisciplinary research area 
in which numbers of researchers with a variety of backgrounds are work- 
ing. An ambition of breaking theoretical boundaries of binary computing 
has been a great driving force behind many investigations to discover quan- 
tum algorithms beyond the boundaries. As celebrated examples beyond the 
boundaries, Shor's algorithms and Grover's one are well-known for the dis- 
crete logarithms, the prime factorization and the data search, respectively 
[H [2]. To trace the history of quantum computing, see for example. 
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On turning to algorithms in convensional sense of computing, there exist 
various excellent algorithms developed in engineering and systems science. 
It might happen that some of them admit a similar mathematical structure: 
In 1990's, Nakamura revealed integrability, Lax-type structure and gradient- 
system structure in a matrix-eigenvalue computing [4j, the gradient system 
on the space of multinomial distributions [S], the Karmarkar flow for linear 
programming [B] and an averaged learning equation of Hebb type . 

In a series of papers by the authors [H [9], [TOl [11] , a counterpart to Naka- 
mura's gradient-system structure was succesfully found in the quantum in- 
formation space (QIS) for the gradient system on the space of multinomial 
distributions |5] and the Karmarkar flow for linear programming [6] . The QIS 
in this paper stands for the space of regular density matrices endowed with 
the quantum SLD-Fisher metric. The former system is realized as the gradi- 
ent system on the QIS associated with the negative von Neumann entropy, 
the latter with a trace of the square of density matrix times a cost-coefficient 
matrix. As a continuation of those papers of the authors, the aim of the 
present paper is to construct a counterpart of the averaged learning equation 
of Hebb type for the principal component analyzer [TJ [T2] , which takes the 
gradient-system form in the QIS. The present paper is placed in an inter- 
disciplinary research area of dynamical systems on the QIS, algorithms and 
applied differential geometry. The result would be expected to be a clue to 
realize algorithms in the QIS in mathematical sense first, and in a more re- 
alistic physical sense in future. In what follows, the contents of the present 
paper are outlined. 

Section 2 is the preliminaries for the averaged learning equation of Hebb 
type (ALEH). The ALEH is derived from a synaptic neuron model and its 
gradient-system form is introduced. Section 3 is devoted to geometric devices 
for realizing the ALEH on the QIS. In order to transfer the ALEH on a dense 
submanifold Sm of the sphere to that on the submanifold T>m of the QIS 
consisting of diagonal matices, an immersion of Sm to Vm is introduced. On 
regarding the immersion as a local diffeomorphism, the ALEH is understood 
to be transferred to a multi-fold copy of the ALHE. It is worth pointing out 
that the geometric devices presented in the present paper are quite different 
from those in the papers [TUl [H] for the Karmarkar fiow. Section 4 is the core 
part of the present paper, where the gradient system on the QIS realizing 
the ALEH is given explicitly. Section 5 is for concluding remarks. 



2 



2 Preliminary: The ALEH 



In this section, we review an averaged learning equation of Hebb-type (ALEH) 
following Nakamura [7] and Oja [12]. The gradient-system form derived by 
Nakamura [7j of the ALEH is also reviewed on the {m — l)-dimensional unit 
sphere. 



2.1 Learning in the synaptic neuron model 

Let us introduce a vector- valued variable X = {Xi,X2, ■ ■ ■ ,Xm)'^ G R™ to 
express the m presynaptic signals and a scalar-valued variable y G R to 
express the postsynaptic signal; the values at the time s of X and Y are 
described as X{s) and Y{s), respectively. The synaptic neuron model dealt 
with in this paper starts with the following timewise linear relation 

Y{s) = W^{s)X{s) (1) 

between X and Y, where W{s) = {Wi{s), W2{s), W„,{s))^ e R™ stands 
for the vector of coupling strengths of the neuron. Throughtout this paper, 
the superscript ^ indicates the transpose operation. 

According to Hebb's hypothesis [121 US], learning in the synaptic neuron 
models amounts to updating efficacies of extracting inputs with high proba- 
bility: In the model ([T]), the vector of coupling strengths W{s) is understood 
to be updated through a recurrence relation along with repetitive inputs of 
signals into neuron. Following Oja [12], we consider the discrete-time recur- 
rence relation 

W{s + l)= ^(^)+vns)x{s) 

^ ^ ^ ||H^(.)+r^y(.)X(.)|| ^ ^ 

for the coupling strength W{s) in ([T]), where rjis a positive constant indicating 
the learning rate and the symbol || ■ || stands for the standard Euclidean norm 
of vectors in R™. As a significant characteristic of the recurrence relation 
([2]), the following norm-preserving property is worth pointed out; 

||H^(.)|| = 1 for ||W^(0)|| = L (3) 

We derive a differential equation providing an approximation of the learn- 
ing system ([1]) with Oja's recurrence rule ([2]) in what follows. Equation ([2]) 
takes the form 

\\Wis) + r]Xis)XT{s)Wis)\\ ^' 
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under the relation ([T]), which admits the Maclaurin expansion 

W{s + 1) = W{s) +r]{X{s)X^{s)W{s) 

- {W''{s)X{s)X''{s)W{s)) W{s)} + 0{rf ) (5) 

if the leaning rate r] is sufficiently small. The 0{rf) in ([5]) denotes the 
second-order infinitesimal. Elimination of the term 0{rf) from the rhs of ([5]) 
therefore provides us with the equation 

W{s + l) = W{s) + r]{X{s)X^{s)W{s) 

-{W^{s)X{s)X^{s)W{s))W{s)}. (6) 

2.2 Averaging 

We proceed an averaging of (jH]) in what follows. Let X{s) and W{s) be 
stochastic processes, which are statistically independent to each other. On 
taking the expectation of ([6]), we obtain 

E[W{s + l)\W{s)]- E[W{s)] 
= v{E[X{s)X''{s)]E[W{s)] 

-{E[W{s)fE[X{s)X^{s)]E[W{s)])E[W{s)]}, (7) 

where the symbol E[-] denotes expectation operation. On assuming the 
stochastic process X{s) to be stationally, the correlation matrix E[X{s)X'^{s)] 
is kept invariant along s, which is thereby diagonalized by an orthogonal ma- 
trix G to 

C = diag(ci, C2, • • ■ , Cm) = G^E[X{s)X^{s)]G. (8) 

Note that the orthogonal matrix G do not depend on s, so that all the 
eigenvalues Cj of E[X{s)X'^{s)] are kept invariant along s, too. The change 
of variables 

w{t) = {wi{t),W2{t), ■ ■ ■ , Wm{t)f = G^W{t/r]) (9) 
with the time-scaling 

t = r]s (10) 

brings ([7]) into the form, 

w{t + r])-w{t) =r](Cw{t) - {w^{t)Gw{t))w{t)y (11) 
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The differential equation 



= Cw- {w^Cw)w (12) 

(JjL 

thereby emerges from fill I) as a continuous-time approximation of (jlj) [3, 112] 
in view of the stochastic approximation theory [13]. Throughout this paper, 
we will refer to f|T2l) as the averaged learning equation of Hebb type (ALEH). 



2.3 The gradient-system form 

We start with showing the norm preserving property 

\\w{t)\\ = l for ||u;(0)|| = l (13) 

of the ALEH (fT2|) . which is understood to be a counterpart to ([3]) of ([1]) with 
Oja's rule ([2]). Indeed, for any solution w{t) of (IT^ with 11^(0)11 = 1, the 
calculation below shows (fT5]l : 

j^\Mt)r = Mtf^it) = 2{w{tfcwmi - Mm = o. (m) 

Putting (P, (UnD and ([HD together, we obtain 

\\W{s)\\ = \\Gw{vs)\\ = \\w{^s)\\ = \\wm = 1 (15) 

for \\W{0)\\ = \\w{0)\\ = 1, as the counterpart to 

Owing to the norm preserving property fll3p . we can restrict the ALEH 
on the (m — 1) -dimensional unit sphere 



w|| = U (16) 



in R"^. We endow 5™ with the standard Riemannian metric 



V {u, u' G T^S'^-\ w G ^"^-^), (17) 



where T^S^ denotes the tangent space of S"*" at w defined to be 



T ,5'"-^ = <i u G R 



w^u = Q\ (we S"^'^] 



According to Nakamura [7], the ALEH on S"^ ^ admits the gradient- 
system form. Let the function, 

A(«;) = -^«;^Ct/7 = -^f;W {weS"'-'), (19) 

k=l 



on 5*™^^ be taken as the potential for the gradient-system form. The gradient 
vector field gradA for the gradient system {S"^~^ , {{■ , , A) is defined as 
follows: For a sufficiently small interval [a, b] with a < < 6, let us associate 
a smooth curve 7 : [a, b] — > S"^~^ with any u G T^jS"^'^ in the manner 



T e [a, b] ^ j{t) e S^^-\ 7(0)=^, ^ 



= ue T^S"'-\ (20) 

r=0 



Then the gradient vector field gradA is defined to satisfy 

A(7(r)) {ueT^S"'-^). (21) 



T = 

Accordingly, gradA turns out to be 

gradA(u;) = -Cw + {w^Cw)w, (22) 



which conicides with the minus of the rhs of f|T2|) . The ALEH on S*™ ^ is 
thus written in the gradient-system form. 

3 Geometric devices 

This section provides geometric devices for realizing the ALEH in the QIS. 
3.1 The QIS 

Following Uwano et al [9], we introduce the quantum information space 
(QIS), the space of regular density matrices endowed with the quantum SLD 
(symmetric logarithmic derivertive) Fisher metric, in what follows. 
Let us consider the space of m x m regular density matrices 

Pm = {p ^ M{m, m)\ = p, trp = 1, p : positive definite}, (23) 

where M{m, m) denotes the set of m x m complex matrices. The Pm is 
endowed with the quantum SLD Fisher metric ((■, ■^'^^ as follows. 
Let the tangent space of Pm at p be defined by 

TpPm = {H G M(m, m) I St = H, trS = 0} . (24) 

The symmetric logarithmic derivertive (SLD) to any tangent vector S G TpPm 
is defined to provide the Hermitean matrix Cp{E?) G M{m,m) subject to 

^{p£,(H) + £,(S)p}=S {EETpPm). (25) 
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The quantum SLD Fisher metric, denoted by ((■, ■))'^^, is then defined to be 
((S, = itr [p{L,{E)L,{E') + L,{E')L,{E))] (S, E' 6 W (26) 

(see also [SI El [HI El [ISl [IS]). 

We wish to present a more ephcit expression of ((■, ■))^^ in what follows. 
Let p G Pm be expressed as 

p = hQh\ h G U(m) 

(27) 

e = diag(^i,...,^„) with trG = 1, 0^ > (A; = 1, 2, ■ ■ -m), 
where U(m) denotes the group of m x m unitary matrices. Expressing S G 

H = /ix/^t (28) 
with /i G U(m) in (1271) . we obtain an explicit expression, 

(h^Cp{:E)h),k = -^^Xjk (j,A; = l,2,---m), (29) 
Uj + Uk 

of the SLD to S G TpPm [9]. Putting (I27D-(I29D into ([26]), we have 

m — / 

{{E,EX^ = 2J2^ (30) 

where E' G Tpi^m is expressed as ([^) 

E' = hx'hl (31) 

The space of m x m regular density matrices, Pm, endowed with the 
quantum SLD Fisher metric ((■, ■))^^ defined above is what we are referring 
to as the quantum information space (QIS) in the present paper, which will 
be denoted also as the pair {Pm, ((■, ■))'^^) henceforth. 

3.2 Metric preserving map 

We start with the Riemannian submanifold 

5^ = ^-^U-iAe\ (32) 

with 

Afj^^ = {w e S'^-'\wk = 0} ik = l,2,--- ,m) (33) 
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of S'"^ ^. From (152]) and ([32]), we immediately obtain the coincidence 

T^S^ = T^S'^-' {weS^cS"'-') (34) 

of the tangent spaces T^Sm and T^S"^~^ (see also ([TS]) ). Then, the metric 
((-, ■))'^ of Sm is defined by 

^0)^ = («, W G T^5„), (35) 

where ^ denotes the differential of the inclusion map 

-.weSm^we 5™-^ (36) 

at w E Sm (see Appendix A for the differential of maps). According to (I5B]) 
and Appendix A, tfy^{u) turns out to be 

'^*,w{u) = ue T^Pm {u G T^Sm)- (37) 
We move to consider the Riemannian submanifold 
Vm = {QePm\Q = diag (^1, ■ ■ ■ , ^™), 

m 

J]efc = l, 4>0 (fc = l,2,--- ,m)} (38) 

A;=l 

of the QIS (Pm, ((■, O)*^^)- The tangent space TeVm at 6 takes the form 

m 

TeVm = [ze M{m, m)\z = diag (Ci, ■ ■ ■ , Cm), J] = o} C T^P^. (39) 

j=i 

The metric ((■, O)^ of is defined by 

((Z, Z'))Z = {{,^^^{Z), i'^^^{Z'))fJ (Z, Z' G Tel)^), (40) 
where denotes the differential of the inclusion map 

: e G ^ e G (41) 

at G Vm (see Appendix A). Like in the case of ^ (see dSZ])), '-^'0(2') turns 
out to be 

l'^^^{Z) = Z eTePm {ZeTeV), (42) 
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Under the preparation above, we consider the smooth map n of Sm to 
of the form 

fi:w = {wi, W2,--- , WmV ^ <Sm ^ diag (wi, wj,--- , u)m) e Dm, (43) 
whose differential of the map /i.,,^^ at w takes the form 

2diag {wiUi, W2U2, ■■■ ,WmUm) {ueTyjSm)- (44) 
We have the following lemma for /i. 

Lemma 3.1. The map /i satifies (i)-(v) in the following: 
(i) The fi is surjective. 

(a) The fi is metric-preserving up to the constant multiple 4; 

4 ((m, u'))^ = ((/i*,^(u), fi'*,w{u')))^(yj) {w e Sm, u, u e T^Sm) (45) 

holds true, where T^Sm is defined by [34\ ) with / fl^) . 
(Hi) Under the {7i2)^-action on Sm defined by 

(f)„ : W = {Wi,W2, - ■ ■ , WmY e Sm 

^ {(yiWi,(J2W2,- ■ ■ ,(JmWmf ^ Sm (46) 

with 

a = ((Ti,a2,---,a„)'^e(Z2)'" (Z2 = {-1,1}), (47) 
the fj, is invariant; 

l2{Mw))=Kw) {weSm,ae{Z2r)- (48) 

(iv) The coincidence ^{w) = ^{w') holds true if and only if there exists a 
certain a G (Z2)™ subject to w' = (j)^{w). 

(v) The restrict of ^ to each of 

S:;^ = {weSm\a,w,>0{j = 1,2,- ■■ ,m)} (a G (Z2)™) (49) 
is difjeomorophic (smooth, infective and surjective) of S"^. 

The items other than (ii) in Lemma 13.11 are proved easily by straightfor- 
ward calculations. The proof of (ii) is consigned in Appendix B. 

Remark 1. If we apply formally to any w G is out of Pm- 

This is the account for eliminating Mrlf^ s from S'^"^ to consider Sm- 

Remark 2. From ( f^^P ([^P with [W) . we have rank/i=K „, = m — 1 = 
dimPm. The map fi is hence called an immersion fTS^ - 

Remark 3. In view of (i) and (iv) in Lemma \3. 1\ T>m is understood as the 
quotient space Sm/ {'^^2)™' of Sm- Further, combining (v) with (i) and (iv), we 
see that iS^ with every a G (Z2)'" is diffeomorphic to Vm- 
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4 The ALEH on the QIS 



We are now in a position to construct the gradient system on the QIS reahzing 
an extention of the ALEH by making full use uses of the geometric devices 
developped in section 3. 



4.1 Mapping the ALEH through fi 



In this subsection, the vector field for the ALEH ([T^ on Sm is shown to be 
mapped to a vector field on Vm- Note that the fact above is not so trivial 
since the map fi is not injective (see (ii) of Lemma [3.11 and Remark [3]). 

In view of (i) of Lemma 13.11 and Remark [H we start with the restrict of 
the ALEH fll2p to Sm- Since gradA satisfies 

(gradA)(u') G T^Af^^ for w e Af^^ (A; = 1, 2, ■ ■ ■ , m), (50) 

all trajectories with initial condition w{0) G Afm^ are confined in 
1,2, ■■■ ,m), respectively. On the same account, all the trajectories with 
w{0) i U'^^-^Min^ never intersect with W^^^Uin^ Accordingly, the ALEH 
m\ can be dealt with on = ^"^-^\ U^^^ Af^\ 

The next thing to show is the (Z2)'"-invariance of gradA, which is equiv- 
alent for 

gradA(0,(w;)) = (0,),,^(gradA(w;)) {w G 5^, a G (Zs)'") (51) 
(see (H6l) and (1471) for 0o-). Indeed, on taking the expression 

{(pa)*,w{u) = ((TiMi, (T2M2, ■ ■ ■ , CTmUm)^ {u G T^^, (J G (Z2)'") (52) 

into account (see Appendix A for the differential of maps). Equation (l5T]) is 
shown to hold true by a simple calculation. 

The (Z2)™-invariance thus shown is put together with (iv) of Lemma [XT] 
to ensure the existence of the vector field /i^fgradA on subject to 

/i^gradA(/i(w)) = ^u^^^ (gradA (w)) {w G Sm)- (53) 

Note that /x*gradA is well-defined since the equation 

/^*,9i.t«(gradA(0^w)) = ;U*,0^^((0<,)^,^(gradA(w;)) 

= /i,,^(gradA(«;)) (ti; G 5^, a G (Z2)™) (54) 

holds true. The first equality follows from fl5T|) and the second one from 

^^*,4>aW O {<Pa)*,w = (Ai O <Pu)*,w = (55) 

with fHHj) . Equation (!55|) follows from fj44|) with fl52|) immediately. 
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Remark 4. On closing this subsection, we give a naive description of the 
vector field /i^,gradA on V^- Recalling Remark we may understand that 
the image T>m of the map fi is a 2"^ -fold copy of Sm- Therefore, in view of 
I^P, the vector field /i^gradA on Dm can be regarded as a 2"^ -folded copy of 
the ALEH on Sm- 



4.2 The gradient system on the QIS reaUzing the ALEH 

We are now in a position to seek a gradient system on the QIS realizing the 
ALEH. Namely, what we we are to seek is a gradient system that realizes 

gradL(e) = /i,gradA(e) (6 G (56) 

where L denotes the potential for the gradient system. We note here that 
the gradient vector field gradL associated with the potential L is defined 
on the QIS (P^, ((-, ■))^-^) in the following way (cf. ^ and ([21])). For a 
sufficiently small interval [a, h] with a < < 6, let us associate a smooth 
curve r : [a, h] — >■ Pm with any S G TpPm in the manner 



T e [a, b] ^ r(r) G Pm, r(0) = p, — 

CLT 



= S G T,Pm. (57) 

r=0 



Then the gradient vector field gradL is defined to satisfy [15] 

L(r(r)) {EeTpPm). (5^ 



(gradL(p),S))?^= A 



T=0 

Instead of (156]) , we are to consider a weaker condition than (156]) below in order 
that we can fix a candidate easily for the potential. The weaker condition to 
be dealt with is 

((gradL(e), if^eZDT = ((/i*gradA(e), Z))g (0 eVm,Ze T^Vm). (59) 
We show the following Lemma. 

Lemma 4.1. // the potential L for a gradient system {Pm, ([■, '))'^^, L) satis- 
fies 

L{fi{w)) = 4A{w) (weSm), (60) 

then l[5U) holds true. 
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Proof It follows from fH5]) that the unique w G S!^ with cr = (1, 1, ■ ■ ■ , 1)-^ 
subject to fi{w) = 9, where is defined by fll^ . Further, recalling (v), we 
can take the unique u G T^Sm subject to = Z. Then for a sufficiently 

small interval [a, 6] with a < < 6, we consider a curve 7(r) subject to (120|) 
and a curve r(t) = /i(7(r)). The setting above is put together with fl2Tl) . 
(SSI), dSSD, dSSD and the assumption (jnUD to show 



(gradL(e),6^e^)) 



e 



(gradL(e),Z)) 



QF 

e 



L(r(r)) 



r=0 



^(/i(7(r))) 



T=0 



4A(7(r)) =4((gradA(«;),M)); 



T=0 



?radA(^i;)),/i,,^(M)))g = ((/i,gradA(e), 



(61) 



This completes the proof. 

Owing to Lemma 14.1^ we can choose 

Lip) = -2tTiCp) 



(62) 



as a candidate for the potential realizing flS^ . where C is the diagonal matrix 
given in (jS]). Note that we can confirm flBU]) for L of fl^ by the calculation, 



L(/i(w)) = -2tr(C/i(w)) 



-2E 



4A(u;), 



(63) 



made with (fT9|) and fH3l) . We move to draw the gradient equation for L of 
(162|) along with the framework of gradient systems on the QIS developped by 
the authors |TT] . Following to pUl [TT] , we define the Hermitean matrix 
;W(L) for L to be 



(A^(L)) 



jk 



dL 

C^jk 

dL 

dp. 



dL 



dpj 



ik 



(l < j < k < m) 

(j = k = l,2,--- ,m). 



where the partial differentiations, d/dpjk and d/dpjk, stand for 



d 



dpjk 
d 



1/ d 



d 



2 \di,u 
If d 



drjjk 
d 



dpjk 2 \d^jk drj 



'jk 



[I < j < k <m) 



(64) 



(65) 
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with 

^jk = ^{pjk), Vjk = '^{pjk) {I < j < k <m). (66) 

The symbols 3? and indicate the real part and the imaginary one, respec- 
tively. In contrast, the pjjS are thought of as real variables to give rise to 
d/dpjjS in usual way. In terms of A1(L), the gradient equation for L is 
written as 

gradL(p) = ^(^pMiL)+M{L)p) - (tr (pA^(L)))p (67) 
(see [ini [H] ) . Since we have 

MiL) = -2C (68) 

through a straightforward calculation of flMl) with fl52]) . the gradient vector 
field for L is exressed as 

gradL(p) = -(pC + Cp) + 2(^tr(pC))p (69) 

with the diagonal matrix C of ([H]). To summarize, we have the following. 

Lemma 4.2. The gradient system [Pm, {{■, ■))'^^ , L) is governed by the dif- 
ferential equation of motion 

^ = (pC + Cp)-2(tr(pC))p, (70) 

where C is the diagonal matrix of 

We are at the final stage to show (1561) for L of (1621) . On combining fl22p 
with (jH]), the rhs of fl56l) is caluculated to be 

p=^,^(gradA(w)) 

= 2diag(wi(gradA(w))i,u;2(gradA(t(;))2, ■ ■ ■ ,ti'm(gradA(w))„) (71) 

with 

m 

(gradA(w))j = -cjWj + ( XI ^kwl^Wj (j = 1, 2, ■ ■ ■ , m), (72) 

fc=i 

while f l69l) is put together with fH3|) to show 

gradL(p(w)) = -2 diag(ciWi, C2wl, ■ ■ ■ , Cmwl^) 

m 

+2 X Cfcw| j diag(ciWi , C2wl, ■ ■ ■ , Cmwl,) (73) 
fc=i 

for w G Sm- Equations fl7T]) - fl73|) therefore confirms fl56|) for L of fl62l) . Finally, 
we reach to the following theorem. 
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Theorem 4.3. The gradient system {Pm, ((■? 'J)^^ iL) governed by the equa- 
tion of motion [71^ realizes the ALEH in the QIS in the following sense: In 
the manner of (E^, the gradient vector field gradL realizes the vector field 
mapped through / f53j) from the gradient vector field describing the ALEH on 

5 Concluding remarks 

We have constructed successfully the gradient system on the QIS (GS-QIS) 
which realizes the ALEH on the submanifold T>m- The success is due to the 
geometric devices developped in Section 3: Especially, a key fact is that the 
immersion /i of Sm to is put together with the (Z2)™-symmetry of the 
ALEH to realize the 2'"-folded copy of the ALEH on V^- 

The behavior of trajectories of the GS-QIS not on V is open still, though 
the trajectories on T) are understood by Nakamura [7j: The GS-QIS is ex- 
pected to have a 'global convergence' property. 

Integrability of the GS-QIS is an open question, too. However, in the 
special case of C = J, integrability is allowed like in the cases of [8] and 
[ini E] with C = —21, since the U(m) action, 

pePm^ hph^ ePm {he U(m)), (74) 

is allowed to be a symmetry, whose role is studied in [8]. 

Relation of the gradient systems obtained in [H [9l [TOl E] to physical 
systems would be a big problem. The horizontal lift of gradient vectors to 
those on the space of ordered-tuples of multi-qubit states developped in P H] 
could be a clue, and the work of Braunstein [16] another one. Further, the 
authors might say that the complementary property between gradient vector 
fields and the Hamiltonian ones in classical mechanics is worth taken together 
with the horizontal lift. 
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oto University for his comments on the present work. This work is partly 
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A Differential of maps 

With the aim of a brief and clear description, we restric to our attention 
to smooth maps among spaces in M(m, m), the space of m x m complex 
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matrices. Let Mi and M2 be manifolds in M[m, m) and a map of Mi to 

M2; 



: Mi^ M2 



(75) 



For a given point p G Mi, the differential V'^^p of ip ai p is defined as follows. 
Like fl20l) . for a sufficiently small interval [a, 6] a < < 6, let us associate a 
sommth curve q : [a,b] Mi 



T E [a, b] ^ q{T) E Mi, g(0) = w, 



dr 



w E T„Mi. 



(76) 



T=0 



On using the smooth curve g(r) introduced above, the differential -^^.^p of the 
map at j9 is defined to be 



^*,p(^) 



d_ 

dr 



^(g(r)) (vETpMi). 



(77) 



r=0 



(see [15], for example). The differential maps, t^*, w, l^q, fi^^w and {(pa)*,w are 
defined by m and dUD with (5„, ^""i, w;, i^), (P^, P^, Q, 6^), (5^, 2^™, w, /x) 
and (5m,5m, w, t/),^) in place of {Mi, M2,p,ip). 



B Proof of Lemma 13.1 



We show (ii) of Lemma 13.11 On setting h = id, S = Z and S' = in ([2 
and (13T|) . Equation (130|) is put together with fHOj) to show 



E 

fc=i 



(78) 



Further the substitutions Z = n^^wiu) and Z' = /i*,u,(M') with (jUj) bring us 
to have 



[Z, z% 



E ¥■ = E = 4 ± = 4((,.. ,/))f.. (79) 



k=l 

This ends the proof. 



fc=i 



k=l 
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